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We search for dimuon decays of a low mass particle in the decays i3° -^ K*°X and 5° — >■ p^X using 
a data sample of 657 x 10'' BB events collected with the Belle detector at the KEKB asymmetric- 
energy e'^e~ collider. We find no evidence for such a particle in the mass range from 212 MeV/c^ to 
300 MeV/c^, and set upper limits on its branching fractions. In particular, we search for a particle 
with a mass of 214.3 MeV/c^ reported by the HyperCP experiment, and obtain upper limits on 
the products B(B° -)■ K*°X) x B{X -^ fJ.+tJ.~) < 2.26 (2.27) x 10"* and B{B° -^ p°X) x B{X -^ 
p'^fj-') < 1.73 (1.73) X 10'* at 90% C.L. for a scalar (vector) X particle. 

PACS numbers: 13.20. He, 12.60.Jv, 12.60.Cn, 12.60.Pr, 14.70. Pw 



The possibility of a weakly interacting light particle 
with a mass from a few MeV to a few GeV has been exten- 
sively discussed [l| . Recent astrophysical observations by 
PAMELA 0] and ATIC [i| have been interpreted as dark 
matter annihilation mediated by a light gauge boson, 
called the C/-boson [J] , which couples to Standard Model 
particles. In addition, the HyperCP collaboration [5[ has 
reported three S"*" — >■ pfj,'^ ^~ events with dimuon invari- 
ant masses clustered around 214.3 MeV/c^ that are con- 
sistent with the process S+ — > pX,X -^ l-i-^l-J-^- Phc- 
nomenologically, X could either be a pseudoscalar or 
an axial-vector particle [6[ with a lifetime for the pseu- 
doscalar case estimated to be about 10"^"* s [7|]. Many 
plausible explanations for such a particle have been pro- 
posed; a pseudoscalar sgoldstino particle [8| in various 
supersymmetric models [9[, a light pseudoscalar Higgs 
boson [lO| in the Next-to-Minimal-Supersyinmetric Stan- 
dard Model as well as a vector t/-boson [ll| as described 
above. 

Recently there have been searches for a similar light 
particle at the Tevatron 
fixed-target experiments 



12 



a 



15| 



^e~ colliders [13[ and 
In those searches, the 
light particle was assumed to be a pseudoscalar and no 
evidence has been found. The KTeV result in Kl de- 
cay disfavors a pseudoscalar explanation of the HyperCP 
results [l5|. 

The large sample of 5° decays at the Belle provides 
a good opportunity to search for a light scalar or vector 
particle. In particular, the estimated branching fractions 
for i3° -^ VX,X -^ A''''M~ where A" is a sgoldstino par- 
ticle with a mass of 214.3 MeV/c^ and V is either a K*'^ 
or p° meson, are in the range 10'^ to 10~^ [16|- 

We report a search for a light particle using the modes, 
B° -^ K*°X, K*° ^ K+TT-, X -)■ u+u- {B%*x) and 



5° -> p°X, pO ^ TT+TT-, X ^ ^+fj,- {B^px) using a data 
sample of 657 x 10^ BB pairs collected with the Belle 



detector [l7| at the KEKB asymmetric-energy e"'"e col- 



lider |l8|. The analysis for B^,^ uses the same dataset 
as Ref. [19|. In this analysis, we assume that the light 
X particle is either a scalar or vector particle. Unless 



specified otherwise, charge-conjugate modes are implied. 
The term scalar (vector) X particle implies either a scalar 
(vector) or pseudoscalar (axial-vector) particle through- 
out this letter. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-flight scintillation counters 
(TOE), and an electromagnetic calorimeter (ECL) com- 
prised of CsI(Tl) crystals located inside a superconduct- 
ing solenoid coil that provides a 1.5 T magnetic field. An 
iron flux-return located outside of the coil is instrumented 
to detect K^ mesons and to identify muons (KLM). 

In the initial event selection, at least two oppositely 
charged muon tracks with momenta larger than 0.690 
GeV/c are required. These muon tracks are selected us- 
ing a likelihood ratio formed from a combination of the 
track penetration depth and hit pattern in the KLM sys- 
tem. We reduce the number of badly reconstructed tracks 
by requiring that \dz\ < 5.0 cm and dr < 1.0 cm, where 
\dz\ and dr are distances of closest approach of a track to 
the interaction point in the beam direction (z) and in the 
transverse plane (r — </)) , respectively. Charged kaons and 
pions are identified using information from the ACC and 
TOE systems and the energy loss {dE/dx) measurements 
in the CDC 0. 

The reconstruction of A'*° {p°) in the B^.^ [B^x) 
decay uses identified A'+ (tt^) and tt' (ti" ) tracks. The 
reconstructed invariant mass Mx*o (Mpo) of K*'^ (p°) 
candidates for the decay mode B^,^ {B^^) is required to 
be in the ranges 0.815 GeV/c^ < Mk-o < 0.975 GeV/c^ 
(0.633 GeV/c2 < MpO < 0.908 GeV/c^), corresponding 
to ±1.5cr (ilcr) in the reconstructed mass distribution. 
The /i^/-t' dimuon tracks are used to reconstruct low 
mass X candidates. 

Br'X (B^x) candidates are reconstructed from a 
j^*o j-pO-j candidate and a pair of muons. Recon- 
structed B^ candidates are selected using the beam- 
energy-constraincd mass ilfbc = \/E^ 
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ergy difference AE ^ Eb - -Bbeam, wfiere £^beam is tlie 
beam energy and Eb {pb) are tfie energy (momentum) of 
the reconstructed B^ candidates evaluated in tlie center- 
of-mass frame. i3° candidates are required to lie in the 
signal regions, 5.27 GeV/c^ < Mbc < 5.29 GeV/c^ and 
-0.03 GeV <l\E< 0.04 GeV (-0.04 GeV < l\E < 0.04 
GeV) for the decay B^,^ {B^-^). In events containing 
more than one B^ candidate, we select the best B^ can- 
didate with the smallest x^ value, where x^ is obtained 
when the four charged tracks are fitted to a common ver- 
tex. Using this algorithm, we select the correct B^,^ 
and B^-^ combinations in the Mbc and AE signal re- 
gion 96.6% (96.7%) and 93.7% (93.5%) of the time for 
a scalar (vector) X particle, respectively. The signature 
for X — )■ /z"*"//^ in i?^*x ^^"^ -^pX decays would be a 
peak in the dimuon mass. The width of the signal region 
for the light particle search with mass below 300 MeV/c^ 
is 3(7 in dimuon mass resolution. The dimuon mass res- 
olutions for B^j^^x and B^-^ vary from 0.5 McV/c^ to 
1.9 MeV/c^ as the mass of X {Mx) increases from 212 
MeV/c^ to 300 MeV/c^. However, the signal region for 
dimuon mass (M^^) for 214.3 MeV/c^ of the HyperCP 
event search is defined to be 211.6 MeV/c^ < M^^ < 
217.2 MeV/c^ where the width of the search region is 
±3(7 in the combined mass resolution, which is obtained 
by linearly summing the mass resolutions of the HypcrCP 
and Belle detectors. 

For background studies, we employ two different tech- 
niques referred to as the counting (C) and fitting {J-) 
methods. Method C uses generic BB and continuum 
(e+e^ ~> qq, q — u, d, s, c) Monte Carlo (MC) samples 
that correspond to an integrated luminosity about three 
times larger than the data sample. In the AE — Mbc sig- 
nal region, there are no events in the dimuon mass region 
M^^ < 225 MeV/c^ (M^^ < 239 MeV/c^) for the decay 
-^K'X i^px)- Ii^ method J", we use the MC samples as 
described above, and select -B" candidates in the side- 
band regions defined as -0.12 GeV < AE < -0.06 GeV 
and 0.06 GeV < AE < 0.12 GeV, and 5.25 GeV/c^ < 
Mbc < 5.27 GeV/c^. By fitting the dimuon mass distri- 
butions for the B^ candidates with a probability density 
function, {x - 0.2iy 



mass in GeV/c 



for X > 2TOp, where a; is a dimuon 
is the niuon mass and the param- 
eter n is extracted from the fit, we estimate the num- 
ber of background events with dimuon mass below 300 
MeV/c^. We also compare the shape of the probability 
density function with the i?° candidates in data sideband 
regions. No significant discrepancy is found. The esti- 
mated numbers of background events for methods C and 
T for the HyperCP event search are (0) and 0.13J;o;q3 
(0.12^0:02) for the decays B^.^ (S°x)' respectively. The 
background estimates for both methods give results that 
are equivalent within statistical errors for masses below 
300 MeV/c2. 

Before examining the full data sample, various distri- 
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FIG. 1: Data and MC comparison for AE and Mhc distribu- 
tions for B^,x (top) and Bpx (bottom) candidates. The 
points with error bars and histograms represent data and 
background MC, respectively. 



butions, including Mbc, AE, dimuon mass and dz in the 
background MC samples arc compared with a small frac- 
tion of the data. These are in good agreement. Figure [T] 
shows the data and MC comparison for AE and Mbc dis- 
tributions after the best B" candidates are selected. The 
peaks in the AE and A/bc distributions for the i?^.x 
are mainly due to B° — >■ J/ijjK*^, J/ij; -^ IJ-'^IJ-~- The 
dimuon mass distributions including the J/ijj and ip' mass 
regions for B^,x and B^x candidates in the signal re- 
gions of Afbc and AE arc shown in Fig. [51 There are no 
events observed in the HypcrCP mass region. 

For the full data sample, no significant signal is ob- 
served for the decays B^,x and B^x for Mx below 
^ 300 MeV/c^. We derive an upper limit for the sig- 
nal yield (S'90) at a 90% confidence level (C.L.) by us- 
ing the POLE program 2l| with the Feldman- Cousins 
method [24I. This procedure takes into account Poisson 
fluctuations in the number of observed signal events and 
Gaussian fluctuations in the estimated number of back- 
ground events as well as systematic uncertainties. The 
Sgo values for the HyperCP event search are 2.33 (2.33) 
for B^,x decay with a scalar (vector) X and 2.33 (2.33) 



for B'^x decay with a scalar (vector) X particle. 

Upper limits on the branching fraction for the decays 
^K'X ^rid B'^x are obtained from 



B{B° ^VX,X ^ H+fi-) < 
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TABLE I: Summary of the number of 
signal yields (Sgo) and upper limits (U. 
errors on Ni,g are statistical only. 



observed events (Nobs), estimated number of background events (Ntg), efficiencies (e), 
L.) at 90% C.L. for the decays B%*x ^nd B'^px with the scalar (vector) X particle. The 
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FIG. 2: Dimuon mass distribution for the B%,x (top) and 
B^px (bottom) candidates in the signal regions for Afbc and 
AU. The shaded region in the inset shows the HyperCP mass 
region. 



where V stands for either K*" or p", and By l23| are the 
intermediate vector meson branching fractions, B{K*^ —^ 
K~^TT~) or i3(p° — ^ 7r"*"7r~). Here Ngg and e denote the 
number of BB pairs and the signal efficiency with small 
data/MC corrections for charged particle identification, 
respectively. 

The signal efficiency is determined by applying the 
same selection criteria to the signal MC sample as those 
used for the data. The signal MC samples for a scalar 
(vector) X particle are generated for X masses in the 



range 212 MeV/c^ < Mx < 300 McV/c^ using the 
P ^VS {P ^ VV) model in the EvtGcn generator ^ 
for a scalar (vector) X particle. In the MC generation 
of the vector X particle, we assume that the polarization 
of X is either fully longitudinal or transverse. The ef- 
ficiency differences between longitudinal and transverse 
polarizations of the X for both modes in the search range 
are less than 7 %. Since the efficiencies for a fully longi- 
tudinal polarized X are lower than for a fully transversly 
polarized X , we conservatively use the efficiencies for full 
longitudinal polarization of the X for upper limit estima- 
tions. In the HyperCP event search for a scalar (vector) 
X particle, the efficiencies for B'j^^x and B^-^ decays are 
23.6% (23.5%) and 20.7% (20.7%), respectively. We also 
check the efficiencies for different X lifetimes. The effi- 
ciencies are the same for lifetimes below 10~^^ s because 
the primary and secondary vertices are indistinguishable. 
The efficiencies for the two different vertex fitting meth- 
ods for the HyperCP event search are compared. One 
method assumes that the dimuon tracks from the X orig- 
inate from the primary S" decay vertex, while the other 
assumes that the dimuon tracks from the X are from 
a secondary vertex. The difference in the efficiencies is 
about 1 %. 

To obtain the final upper limit, we use the backgrounds 
determined from the fitting method. Since the efficien- 
cies for a scalar (vector) and a pseudoscalar (axial-vector) 
are the same, the upper limits for the scalar (vector) 
and the pseudoscalar (axial-vector) X searches are iden- 
tical. From the B^,^ (-^px) sample, the upper lim- 
its for a scalar and vector X particle in the HyperCP 
mass range are determined to be 2.26 (1-73) x 10~^ and 
2.27 (1.73) X 10~^, respectively. Table |T] summarizes the 
number of observed events, the expected number of back- 
ground events, the efficiencies, the signal yields, and the 
upper hmits at 90% C.L. in the interval 212 MeV/c^ < 



TABLE II: Summary of fractional systematic uncertainties in 
the upper limit for a scalar (vector) X particle in the HyperCP 
mass range for the decays -B^.^ and B^xj respectively. 



Source 


ae/B (%) 


rO 
Bpx 


[IJ 


Nbb 


1.4 (1.4) 


1.4 (1.4) 




fi^ identification 


4.2 (4.2) 


4.1 (4.1) 




K"^ identification 


0.8 (0.8) 


- 




n identification 


0.5 (0.5) 


1.0 (1.0) 


[A 


Tracking efficiency 


4.2 (4.2) 


4.4 (4.3) 


[3] 


Afbc 


0.5 (0.3) 


0.3 (0.6) 


AE 


0.5 (0.3) 


0.3 (0.6) 


[41 


K*° tagging 


0.5 (0.3) 


- 




p° tagging 


- 


0.3 (0.6) 




MC statistics 


0.1 (0.1) 


0.1 (0.1) 


[5] 


Total 


6.2 (6.2) 


6.2 (6.3) 


161 



Mx < 300 McV/c2. 

The systematic uncertainties in the upper limits for the 
decays i?/^.^ ^^^ ^'pX i^^ ^^^ HyperCP mass range are 
summarized in Table |lTl The total systematic uncertain- 
ties in the upper limits for both decay modes vary from 
6% to 8% as the mass of X increases from 212 MeV/c^ 
to 300 MeV/c^. The dominant systematic uncertainties 
come from tracking efficiency and muon identification. 
The uncertainty for the tracking efficiency is estimated 
by linearly summing the single track systematic errors, 
which are ~ 1%/track. The uncertainty of muon iden- 
tification is measured as a function of momentum and 
direction by using the 77 — >■ /i^/i~ data sample. 

In summary, we searched for a scalar and vector par- 
ticle in the decays B° -^ K*^X, K*° ^ A'+tt", X -^ 
/^+/x~ and i?° — )► p°X, p" — > tt+tt^, X — > fJ,^ p^ in the 
mass region 212 MeV/c^ < Mx < 300 MeV/c^. No sig- 
nificant signals are observed in a sample of 657 x 10^ BB 
pairs. We set 90% C.L. upper limits of B{B° -^ 
K*°X, K*° -^ K+TT-, X -^ p+p^) < 2.26 x lO"® 
(2.27 X 10-8) and BiB° -^ p°X, p° -^ tt+tt", X -^ 
p+IJL-) < 1.73 X 10-8 (1.73 X 10-8) for a 214.3 MeV/c^ 
mass scalar (vector) X particle; our results rule out mod- 
els II and III for the sgoldstino interpretation of the Hy- 
perCP observation [16|. 
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